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Abstract: We study the effect of departures from SUSY GUT universality on tlie 
neutralino relic density, and both its direct and indirect detection, especially by 
neutrino telescopes. We find that the most interesting models are those with a value 
of M^Igut lower than the universal case. 
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1. Introduction - CMSSM summary 

In a supersymmetric framework with i?— parity conservation, the lightest supersym- 
metric particle (LSP) is stable. In the Minimal Supersymmetric Standard Model 
(MSSM), it is often the lightest neutralino (= the neutralino x) which is a neu- 
tral Majorana particle. It then offers an interesting candidate to account for cold 
dark matter (CDM) in the present Universe {^cdm ~ 0.3). The relic population of 
these neutralinos which survive after spatial separation freezes their self-annihilation, 
could be detected by the energy they transfer to nuclei in direct detection experi- 
ments. Another possibility is the indirect detection of fluxes coming from the decays 
of neutralino annihilation products {xx ~^ ~^ ^ili^^ •••) (fo^' ^ review on neu- 
tralino dark matter and different detection possibilities see [0]). Since these fluxes go 
like the squared neutralino density, they require some concentration to restart their 
annihilation. The gravitational concentration achieved in galactic halos is much too 
low for neutrino indirect detection. However, neutralinos can also accumulate in 
the gravitational well of massive astrophysical bodies, and the larger concentration 
achieved in the centre of the Sun can give rise to higher muon neutrino fluxes de- 
tectable by neutrino telescopes such as Antares or IceCube 0. 

In a previous paper |Q we studied the potential detection of neutralino dark mat- 
ter by neutrino telescopes and direct detection experiments in CMSSM models (also 
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known as mSugra). Those models inspired from minimal supergravity and gravity- 
mediated SUSY breaking assume a unification of the soft parameters of the MSSM 
at high energy Mqut ~ 2 x 10^^ GeV reducing the 106 "SUSY" parameters of the 
MSSM down to 5 : universal masses for scalars (mo) and gauginos (mi/2), univer- 
sal trilinear (Aq) and bilinear (i?o) couplings, and a Higgs "mass" parameter (yUo)- 
Using renormalisation group equations (RGE) and requiring radiative electroweak 
supersymmetry breaking, the usual input parameters are mo, mi/2, Aq, tan/? (the 
ratio of the 2 Higgs doublet vacuum expected values at low energy) and sign(/i). 
The neutralinos are the mass eigenstates coming from the mixing of neutral gauge 
and Higgs boson superpartners. In the MSSM the neutralino mass matrix in the 
{B,W\H^,Hl) basis is : 



^ Ml —mz cos (3 sin 9^r mz sin /? sin 6^r ^ 

M2 mz cos P cos —mz sin {3 cos 6^r 

^ —mz cos [3 sin 6-^ mz cos [3 cos -fi 

y mz sin /5 sin 6y^ —mz sin [3 cos 6-^ —fi J 

(1.1) 

and can be diagonalised by a single mixing matrix z : Mdiag = zMnz~^. The 
(lightest) neutralino is then given by the linear combination 

X = ZnB + ZuW^ + z,,H^, + ^14^2°- (1-2) 

In CMSSM models, the neutralino exhibits two different natures: it is fully bino- 
like (zii ~ 1) for low mo values because of RGE evolution of parameters in eq. ( |1.1|) , 
and can acquire a non-negligible higgsino component for mo > 1000 GeV along the 
boundary where radiative electroweak symmetry breaking cannot be achieved be- 
cause of the focus point behaviour In a previous work f^, we used Suspect^ [|, ^ 
for RGE, potential minimisation and SUSY spectrum calculations, and DarkSusy[|] 
to estimate relic density and detection rates. We carefully analysed the dominant 
annihilation channels which count both in relic density calculation and neutrino 
spectra. In regions where the neutralino relic density satisfy current cosmological 
constraints, we found that: 

• the low mo and mi/2 region is strongly constrained by experimental limits on 
SUSY contributions to the 6^5 + 7 branching ratio, on the lightest Higgs 
mass and on SUSY contributions to {g — 2)^ : a^^^"^ , 

• xt coannihilation regions are beyond reach of detection, 

• the pseudo-scalar A pole region can only be potentially detected with very big 
future projects for direct detection (~ 1 ton size). 



^Suspect has been recently updated. This improvement does not affect our results. 
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mSugra with 5 GeV threshold vs neutrino Indirect Dection 



mSugra models vs Direct Detection 
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Figure 1: a) Neutrino indirect detection experimental (Macro [Q], Baksan [10|, Super- 
Kamiokande Antares |12], IceCube sensitivities on muon fluxes with a 5 GeV 
threshold coming from x annihilations in the Sun and b) direct detection experimental 
(CDMS HI, Edelweiss l|5| and II [||], Zephn Q) sensitivities on cr^™^. Both are function 
of the neutralino mass for a wide scan of CMSSM models (see text for details). 



• the large mo "focus point" region where the neutralino has a significant (and 
crucial, see figure ) higgsino fraction (///) is interesting for direct and indirect 
detection experiments. Neutrino/muon fluxes coming from the Sun are large 

only in this region because of the XX — ^ W~^W~, ZZ and XX — ^ channels 
which give rise to more energetic neutrino spectra and to muons with higher 
energy. 

This is summarised on figures |I| and ^ for the following sample of CMSSM models, 
compared to sensitivities of current and future direct and neutrino indirect detection 
experiments : 

• < mo < 3000 GeV ; 40 < mi/2 < 800 GeV ; Aq = GeV ; tan /3 = 10, 50 ; 
/i > 0, 

• < mo < 3000 GeV ; 40 < my2 < 1000 GeV ; Aq = -800, -400, 0, 400, 800 
GeV ; tan/? = 20, 35 ; > 0, 

• < mo < 3000 GeV ; 40 < mi/2 < 1000 GeV ; Ao = GeV ; tan /5 = 45 ; 
/X > 0. 

The constraints we apply on the models are : 
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Figure 2: a) Neutralino relic density and b) muon flux coming from the Sun as functions 
of the X higgsino fraction for a wide scan of CMSSM models (see text for details). The 
band in plot a) indicates the current preferred range for VLcdm- 



limits on SUSY particle masses from accelerators searches (see e.g. |18| to 
update these values): 

> 104 GeV; m^- > 100 GeV for / = ii,bij^,9, rrig > 300 GeV; m,-, ^ > 
260 GeV for q = u,d, s, c, 

limits on Higgs mass [|1^] : we require rrih > 112 GeV instead of m/j > 114 
GeV (depending on tan(/3 — a)) because of the too low value for Higgs mass in 
Suspect 2.005 (~ 3 GeV) as pointed out in |20[| , 



limits on 6 ^ S7 branching ratio : we require BR{b ^s + 7) = 1^4x 10~^ 
as calculated in DarkSusy, corresponding to BR{h s + 7) ~ 2.3 5.3 x 10~^ 
(the Standard Model value calculated in DarkSusy is 2.4 x 10~^ instead of the 
NLO 3.6 X 10~^ 1^) when experimental results are BR{b — s + 7) = 
3.37 ± 0.37 ± 0.34 ± 0.24 ± 0.38 x 10"^ 



• limits on SUSY contribution to muon anomalous moment [g — 2)^ |^3| : —6 x 

10-10 < aff^^^ < 58 X 10-10. 

The Renormalisation Group Equations (RGE) lead SUSY models to a generic 
hierarchy of particle spectrum in which scalars are heavier than light neutralinos 
and charginos. As summarised above, the most interesting CMSS Models for the 
detection of a relic neutralino are those for which its higgsino fraction is non negligible 
and has a dominant effect (figure These are indeed the only models leading to 
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large neutralino annihilation cross section (important for its relic density and 
for indirect detection) and to large neutralino-proton scalar and spin dependent 
elastic cross sections (important for direct and indirect detection). In this paper, 
using the same tools as previously f^, |[ , we shall relax some universality hypotheses 
and examine whether more favourable models for relic density and detections (both 
direct and neutrino telescopes) can be found, and establish that models detectable 
by neutrino telescopes are more generic than the "focus point" region of a typical 
(mo,mi/2) CMSSM plane, offering a less constrained framework for detection. We 
will then conclude by a generic low energy parameterisation of typical models coming 
from RGE evolutions, favourable for neutrino indirect detection. 



2. Non-universal SUSY GUT models 

In which direction can one relax the universality hypothesis in order to induce the 
physics we want, namely to get large cr^'^g and c^lq and thus a non negligible 

higgsino fraction for the neutralino ? 

We can start with some arguments on non-universality inspired by the (1-loop) 
renormalisation group equations p4|. It has been shown pSf that from the RGE 



solutions of the couplings in the unbroken symmetry phase, one can obtain the soft 
terms of the broken phase by an expansion over Grassman variables. This has been 
done for the MSSM couplings ||2^ in order to get the RGE solutions for the soft 
terms. In the case of an analysis with three independent Yukawa couplings ht, hi, 
and hr, it is found that due to the infra-red quasi fixed point (IRQFP) [|2^ behaviour 
of the Yukawa couplings, the low energy values of the third generation scalar soft 
masses and of the Higgs masses have a weak (scale independent) relation with their 
initial values at Mqut and depend mainly on the high energy 5*^/(3) gaugino soft 
mass M^\gut (in a scale dependent way). The first and second generation soft masses 
have an analogous behaviour because of their negligible Yukawa couphngs. This can 
be written as : 

{MZjWio^f = {MZJWgut? + C3/3 + C2/2 + ci/i + corrections (2.1) 

with 

h \ (1 + W)V' * ' 

where Mj are the soft gaugino mass terms and is the universal gaugino coupling at 
Mguti and where C3 is strongly dominant for squarks and Higgses and less dominant 
for sleptons which have vanishing SU{?)) charges. 

In addition, a gaugino non- universality given by {M2/ Mi)gut < 1 can lead to 
a large wino component for the neutralino and to an important modification of its 
couplings with respect to the universal CMSSM case and thus to a very different 
phenomenology. 
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We will now explore various non-universal scenarios in order to study their pos- 
sible benefits on detection. 



3. Non-universality of scalar soft masses and of trilinear cou- 
plings at McuT 

Sfermions : 

As for Aq in the CMSSM modifications of the soft parameters in the sfermions 
mass matrices can give rise to light third generation sfermions. This could modify the 
neutralino relic density value through coannihilation processes. However, in detection 
processes, the neutralino interacts with nuclei and thus mainly with u and d valence 
quarks. Due to their small Yukawa couplings, RGE evolutions of the first and second 
generation squark masses only depend on gaugino soft masses. This implies that their 
masses can not be lowered by changing scalar soft terms to enhance cr^'^^ and cr^l^ 

through the process XQ XQ- Non-universality in the sfermion mass soft terms can 



thus only lead to ' Xt PH) coannihilation effects, giving rise to models 

for which the neutralino is cosmologically a good dark matter candidate without 
modifying its detectability. 

The soft terms for the third generation sfermions being present in soft Higgs 
masses RGE's (Xj , Xf, and X^) can also influence the radiative electroweak symmetry 
breaking through the running of m|^^ and m^H2- This effect is similar to the one 
discussed in the next section. One could also add lighter third generation squarks 
but this is not helpful for detection. 
Higgses : 

One can relax the universality relation for the Higgs soft masses at the Grand Uni- 
fication scale mHilcuT = i^h-zIgut = mo- and mj^j being now free parameters, 
the potential minimisation and thus the value of /i is less constrained. The parame- 



terisation usually taken is ||3T|, |32|, |33 



m h,\gut = {I + Si)mo ; for i = 1,2. (3.1) 

Due to the departure of the Higgs soft masses with respect to the universal case, /i 
values are modified and the parameter space of a typical (ttiq, ""^1/2) plane satisfying 
the radiative electroweak symmetry breaking is reduced with respect to the CMSSM. 
In the remaining parameter space, models with smaller masses for the heavier Higgs 
bosons TTiA and tuh are more easily obtained than for the (universal) CMSSM, lead- 
ing to wider zones with a good neutralino relic density and accessible direct detection 
yields. Concerning the neutrino indirect detection potential, the muon fluxes due to 
neutralino annihilation in the Sun remain however small for such models. This is 
in contradiction with the recent results presented in ref. [Q, in which the modifica- 
tion of the Higgs soft mass relations generate more easily models with large muon 



6 



= niQ ', m-Hi = 0.5mo ; Aq = ; tan /? = 45 ; jj, > 




500 1000 1500 2000 2500 3000 500 1000 1500 2000 2500 3000 

mO (GeV) mO (GeV) 



a) b) 

Figure 3: Neutralino detection potential a) on muon fluxes coming from the Sun for 
neutrino telescopes and b) for direct detection experiments in the (jti-Oj "1-1/2) plane for 
models with non-universal Higgs soft masses rn'jj^ = niQ and rnj^^ = 0.5mo and with ^0 = 0, 
tan/3 = 45, > 0. Also shown are constant relic density lines ^y^h"^ = 1, 0.3, 0.1, 0.025. 
The small grey area shows the models excluded by current experimental constraints. 



fluxes coming from the Sun in the parameter space allowed by radiative electroweak 
symmetry breaking. This difference can certainly be attributed to the use in 
of IsaSusy (version < 7.64) for RGE's and SUSY spectrum calculation (instead of 
Suspect in this work) for which the predicted value of fi is usually too small for high 
mo values |20|, and hence the higgsino fraction of the neutralino is too large. An 
example of our results is shown on figure ^for 61 = —0.5 and 62 = 0. 



4. Non-universality of gaugino soft masses at Mqut 
4.1 SU{h) representations 

In CMSS Models, a universal mass is given to all gaugino fields at the high energy 
scale Mgut- The unification gauge group should either be SU{^) or a larger group 
naturally breaking into an SU{b) subgroup at a scale around Mqut- In the class of 
models for which the SUSY breaking is induced by an F-term, the gaugino masses 
are generated by a chiral superfield with an auxiliary component acquiring a vev 
such that H, |36|, ^ 

Mgauginos ~ —X"'X'' (4.1) 

^y-lPlanck 
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Table 1: Relative values of the gaugino masses at Mqut and mz scales in the four possible 
irreducible representations for F^. 



where A" and are the gaugino fields B, W and g. Since all gaugino fields belong 
to the adjoint representation, F$ belongs to an irreducible representation in the 
symmetrised product of two adjoints of SU (5) (=24) or to a linear combination of 
these representations : 

(24 X 24)sym = 1 © 24 © 75 © 200 (4.2) 

where only the singlet component (1) leads to universal masses Mi\gut = M2\gut = 
MsIgut = mi/2. The embedding coefficients of the Standard Model gauge groups in 
SU{5) give the relations between the gaugino masses at Mgut (see table [^), which 
we now discuss for the pure non-singlet irreducible representations. 
Case of the 24 : 

For low mo values, since MfV^I^ < M^^^^^ /M§^^^^^ ~ 1/2 at low energy the 
neutralino is strongly bino-like. As for the CMSSM, this kind of models will not 
lead to large indirect detection muon fiuxes coming from the Sun. In addition, the 
lower limits on the neutralino mass obtained by LEP experiments are only valid for 
M1/M2 > 1/2 at low energy, so models with a very light neutralino are not excluded. 
But the main grey bubble part of the exclusion domain comes from the lower limits 
on the SUSY contribution to the muon g — 2 parameter a^^^ which are very strong 
for these models. 

For large mo values, a small higgsino fraction in the neutralino can be obtained 
in the focus point region, but less easily than in the CMSSM universal case because 
the larger IM2I value leads to larger m^^^ values. So the direct detection yields 
and indirect detection muon fiuxes for these models are less important than for the 
universal case and the neutralino relic density is too large. This is illustrated in 
figure where no model is accessible to the Antares sensitivity, and where the very 
thin IceCube detection region always leads to a cosmic over- weight. In addition, the 
Edelweiss II and Zephn detection regions are strongly reduced with respect to the 
universal case for large mo values, due to the small higgsino fraction of the neutralino, 
and also due to the relative sign between Mi, M2 and (same effect as yU < for 
the CMSSM with the "hole" in the spin independent elastic cross section due to the 



suppression between up and down quarks contributions |^ 
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Figure 4: Neutralino detection potential a) for neutrino telescopes detecting muon fluxes 
coming from the Sun and b) for direct detection experiments in the (mo,— Mi) plane of 
the 24 representation for = 0, tan/3 = 45, /x > 0. No model with a good relic density 
is accessible to Antares (0.1 km^). 



This is in agreement with the results of ref. in which the neutrahno has been 



chosen even more bino-hke, M2/M1 ~ 10, leading to a relic over-density except for 
light sleptons. 
Case of the 75 : 

At low energy \Mi\ ~ IM2I ~ IM3I, the neutralino is thus equally wino and bino, 
and the two lightest neutralinos and the lightest chargino are almost degenerate 

~ M^o ~ M^+. These models open up the very efficient XX ~ — ^ W~^W~, ZZ 
annihilation channels, as well as very strong XX2 ^-nd xxt coannihilations. This leads 
to a strong reduction of the relic density. In addition, the gluino can sometimes be 
the LSP. 

Case of the 200 : 

Here IM2I ~ 2/5|Mi| at low energy giving a strongly wino-like neutralino. xX ~~ — ^ 
W^W^, ZZ annihilations and xxt coannihilation totally suppress the neutralino 
population. 

4.2 Free relations in gaugino mass parameters : effect of M^Iqut 

As explained above, the MSSM non-universality parameters giving the strongest 
impact on the detection of dark matter neutralinos are M2\gut and mostly M^Igut- 
We will now study the departure from universality of these two parameters and their 
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benefits on the neutralino relic density and the detection yields. These non-universal 
values will then be translated in the above SU{5) representation decomposition. 



mo = 1500 ; mi/2 = 600 ; Aq = ; tan/3 = 45 ; /x > 
a) b) c) 




d) e) f) 

Figure 5: Evolution of a) the neutralino relic density, b) the muon flux coming from the 
Sun, c) the spin independent neutralino-proton cross section (direct detection), d) the fj, 
parameter, e) the gaugino fraction and f) the pseudo-scalar mass ttia (f) as functions of 
the and ^ ratios for the CMSS Model with mo = 1500 GeV, mi/2 = 600 GeV, 
Ao = GeV, tan/3 = 45, /i > 0. 



In the following, the departure from universality effects are quantified by the 
ratios at Mqut x = M2/ 1711/2 or x = M^/mi/2 (with mi/2 = Mi = M3 or mi/2 = 
Ml = M2 respectively) which will be lowered starting from the CMSSM case {x = 1). 
The M2\gut parameter: 

The effect of the M2 parameter is essentially a modification of the neutralino com- 

position. When the wino component of the neutralino increases, the XX — — ^ 
W~^W~ , ZZ processes become more effective and enhance the annihilation cross 



section cr^_^ [0. In addition, the strong XX2 Xxt coannihilations become ac- 
tive and the neutralino relic density strongly decreases. This wino component can 
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increase detection rates by an order of magnitude at most. Indeed, the neutralino- 
quark coupling, which enters in direct detection {(J^^q, H exchange) and in the 
capture for indirect detection {cr^lq, q exchange), is tan 6'iy-suppressed for pure bino 
w.r.t. pure wino. The neutrahno annihilations into the hard W^W~ spectrum also 
give rise to more energetic muons. These enhancements of the (in) direct detection 
yields can reach several order of magnitude with respect to the CMSSM case depend- 
ing on the values of the other SUSY parameters as can be seen in figures ^ and 
^ and l^c, and |^ and |^. 

However, the relevant value of M2 is very critical, so its benefits are only oper- 
ative in a very narrow range. Given an M2/mi/2 ~ 0.6 — 0.7 ratio (equivalent to 
Mi\iow ~ M2\iow)^ the neutralino detection yields are enhanced but the relic density 
drops down to much too small values (see figures |^, |^a and ^ ). In conclusion, 
the handling of this M2 parameter in order to get the desired neutralino dark mat- 
ter phenomenology can only be done by "fine-tuning". One way around this wino 
neutralino extermination suggested in is to have the wino neutralino population 
derived from AMSB models {Mi\iow — 3M2|/o«)) regenerated at low temperature by 
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moduli decays which could give a good relic abundance. 
The MsIgi/t parameter: 

The impact of variations in the M3 parameter is much more interesting. It is indeed 
one of the key parameters of the MSSM through the RGE's. Its influence goes well 
beyond the neutralino sector. Indeed, following the RGE (see equation p?T|) , 
a decrease of M^\gut leads to a decrease of m|^^ and of /i through the radiative 
electroweak symmetry breaking mechanism, thus enhancing the neutralino higgsino 
fraction, and leads also to a decrease of iriq and m^. These effects are illustrated on 
the figures ^J^,^, ^d,|^J^f and ^J^,^. The neutralino relic density then gradually 
decreases with x = M3/mi/2 (see figures |^a, |a and ^ ) due to the increase of the 
CMSSM dominant annihilation cross section channel (mainly XX ~* bb, XX — ^ 

XX W^W~ and xX ZZ according to the CMSSM starting parameters |^). 
The annihilation channels which directly depend on the neutralino higgsino fraction 
finally dominate when x is further lowered, because of the decreasing of fi. 

For the CMSS Model with mo = 1500 GeV, mi/2 = 600 GeV, Aq = GeV, 
tan P = 45, > (figures ^ and ||a ), the dominant channel at x = 1 is xx By 
decreasing x, the latter remains at first dominant while the neutralino relic density is 
reduced due to the decrease of (figure ^ and f ), then since also decreases, the 

processes xX ~^ followed by xX W^W~ and xX take successively over 

(as well as the XX^ ^ind XX2 coannihilations) further lowering the relic density. The 
enhancement of the tt branching ratio before the increase of the neutralino higgsino 
fraction (x > 0.8) (see figure ||a ) is due to the decrease of the stop mass. 

For the same CMSS Model but with tan/5 = 10 (figures |^ and ^ ), is larger 
due to the smaller tan/3 value, so the xX — ^ channel dominates the neutralino 
annihilation. When x is lowered, this process remains dominant (figure ) but its 
cross section firstly increases because decreases. Then the Z exchange process 
proportional to the neutralino higgsino fraction takes over due to the decrease of 
/i (see figure ^ ) enhancing the annihilation cross section and decreasing the relic 
density (see figure |]a ). Finally, the annihilation processes into gauge bosons become 
dominant for x < 0.5 because becomes smaller than the top mass. 

For the CMSS Model with mo = 3000 GeV, mi/2 = 2000 GeV, Ao = GeV, 
tan/3 = 45, /i > (figures |^ and ^c ), the neutralino relic density is much too large 
in the CMSSM, but lowering M3 drives this model into the cosmologically favoured 
region by the decrease of and fi (figures |^, 0d and |^ ). This model is located 

further away from the radiative electroweak symmetry breaking boundary and has 

— 

larger /i. In addition, < m,q so the process xX ~^ strongly dominates over 

z t — 

XX tt. Again, annihilations into gauge bosons take over when x {i.e. /i) decreases. 
Neutralino direct and indirect detections: 

When X is lowered, the direct detection yields (see figures ||c, |^c and |^ ) can be 
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Figure 7 : Same as figure ^ for tlie extreme CMSS Model with mo = 3000 GeV, 
2000 GeV, Ao = GeV, tan/3 = 45, /i > 0. 
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Figure 8: Evolution of the dominant neutralino annihilation branching ratios as functions 
of the Ms/mi/2 ratio for the models a) of figures |5|, b) of figures ^ and c) of figures 0. 



13 



enhanced by 2 or 3 orders of magnitude with respect to the CMSSM case {x = 1). 
Firstly the reduction of the squark masses favours the XQ XQ process, and mainly 
the coupling C^qH is maximal for maximum ^i3(4)-2ii{2) mixed products which increase 
with the neutralino higgsino fraction when x decreases. Moreover, ttlh ~ ttla and 
thus cr^™p is enhanced due to the decrease of m^. When the neutralino gaugino 
fraction finally drops, the 2:13(4)^11(2) products decrease and decreases back with 
X. This behaviour can clearly be remarked on the figures |^, |^, ||e, ||f and 

0c, 0e, ^. 

As far as neutrino indirect detection is concerned, the enhancement on the muon 
fluxes coming from neutralino annihilation in the Sun, due to the decrease of M^Igut, 
can reach up to 6 orders of magnitude with respect to the CMSSM case (x = 1). 
The main effect is coming from the increase of the spin dependent neutralino-proton 
elastic cross section, firstly due to the decrease of the squark masses which enhance 
Xq XI ill '^x-p' then mainly to the decrease of leading to a larger higgsino 
fraction in the neutralino which enhance XQ ~* XQ '^x-p- Moreover the larger hig- 
gsino fraction also favours the neutralino annihilations into the XX ~^ W~^W~ , ZZ 
and XX ~^ tt channels which give harder neutrino spectra than xX ~^ bb. This is 
illustrated on figures ^b, ||d, ^e, ||b, ^e and 0b, 0d, 0e. This enhancement is not 
as peaked in x as for direct detection but remains maximum as long as the higgsino 
fraction dominates the neutralino composition. However the relic density becomes 
very small when x is further lowered. For these extreme low values of M3 with small 
relic density, the enhancement of the muon fluxes coming from neutralino annihila- 
tion in the centre of the Earth (enhanced by cr^^p and o"^.^) can become significant 
but usually remains beyond reach of current and next generation neutrino telescopes' 
sensitivities. 

To summarise, we have shown it is possible, by lowering M^Igut, to decrease the 
neutralino relic density to the desired cosmological value for any CMSS Model even 
for mo and mi/2 as large as several TeV. The value of M^Igut necessary to get a relic 
density f^^^^ ~ ^.1 — 0.2 mainly depends on mi/2 : 0.5mi/2 < M3 < mi/2- 

In the CMSSM, the main neutralino annihilation channels can be regrouped in 

— 

two sets: the pseudo-scalar exchange xX ~^ bb and the processes directly proportional 

to the neutralino higgsino fraction xX ~^ XX W^W~ and xX ZZ. As we 
have seen, the decrease of the M^\gut parameter influences many MSSM parameters 
through the RGE and leads to both a decrease of and an increase of the neutralino 
higgsino fraction by the decrease of /x. So the dominant neutralino annihilation 
process of any CMSS Model can be enhanced in order to obtain the good value for 
the relic density. 

Varying the x value allows to have f^^/i^ ~ 0.15 in the whole shown (mo,mi/2) 
plane. For tan (3 < 50, where the A pole is not present in the neutralino annihilation, 
the interesting values of x range between ~ 0.8 — 1 along the zone with a correct 
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Ao = ; tan /? = 45 ; /i > ; Mg/mi/s = 0.63 




i:a.^ipie.^i;a^g^. I I I I I I I I I I I I I I I I I III!,,,, i^p^^^^^^^^^^^^^^^ 
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a) b) 

Figure 9: Neutralino detection potential in the ("i0)"^i/2) plane for x = M3/mi/2 = 
0.63 for a) neutrino telescopes looking at muon fluxes coming from the Sun and b) direct 
detection experiments. 



neutralino relic density in the CMSSM and down to x ~ 0.5 — 0.6 for large mi/2 > 
2000 GeV and mo < 2000 GeV. At the centre of the shown (mo,mi/2) plane, a quite 
generic x value is found 

M^Igut ~ 0.6(±0.1)mi/2 + corrections(mo, tan/5, mft) (4.3) 

permitting to get close to the A pole or at least to adjust the couple (m^, -2ii(2)-2i3(4)) 

A — 

and thus the xx — ^ bb channel in order to get ~ 0.1 — 0.2. These x values 

also favour the direct detection yields through the XQ — ^ XQ process. Moreover, 
the decrease of M3 increases the neutralino higgsino fraction, so that a typical value 
^3/^1/2 ^ 0.8 enlarges the "focus-point" corridor at large mo as well as the region 
of the parameter space accessible to direct and indirect detection. Notice that all x 
values discussed above lie well above x = 0.16, below which the gluino becomes the 
LSP. 

Two examples with a fixed ratio M^/mii2 are shown on figures |^ and |10| where 
regions with interesting relic density and experiment sensitivity areas are vastly im- 
proved with respect to the CMSSM case [Q, especially noticing the larger mi/2 range 
shown. 

Links to SUSY breaking : 

From the above discussions, M^\gut appears as the most relevant parameter in SUSY 
non-universal models as far as neutralino detection is concerned and we have defined 
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Ao = ; tan /? = 10 ; /i > ; Mg/mi/s = 0.55 
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Figure 10: Neutralino detection potential in the {mQ,mi/2) plane for tan/3 = 10 with 
X = M^/mi/2 = 0.55 for a) neutrino telescopes looking at muon fluxes coming from the 
Sun and b) direct detection experiments (For this value of x, Zeplin Max covers the whole 
shown plane). 



its departure from universality by the x = gut /''t^'I/ 2 ratio. In line with the logics 
of minimising the deviations to universality, we can map this single parameter x into 
a single combination of the 24, 75 and 200 representations: 

non-univ = CMSSM + C2424 + 075X5 + C2oo200 (4.4) 

which leads for the relative values of the gaugino masses (Mi : M2 : M3) at Mqut to 

(1:1: x)non univ = (1:1: l)cMSSM + (0 : : a; — l)c24 24+C75 75+C2oo200- (4.5) 

From table [|, we then see that the coefficients must take the values: 

C24 = ^(x - 1) ; C75 = ^{x - 1) ; C200 = ^{x - 1). (4.6) 
63 63 63 

giving for the typical value of interest x = 0.6 : C24 — —0.13 ; C75 ~ —0.09 ; 
C200 — —0.06, corresponding to quite small coefficients indeed. 

A possible origin of such lower M3 values might otherwise be found in anomaly 
mediated SUSY breaking (AMSB) schemes. Even if the usual derivation of AMSB 



spectra leads to wino LSP's 0, other patterns have been found where the 
gluino can become lighter, to the point of having a gluino LSP in extreme cases. An 
exploration of this connection in the general framework of MM would be interesting. 
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5. Low energy effective parameterisation of the MSSM 



We have explored the grand unification SUSY models favourable to the indirect 
detection of neutralino dark matter with neutrino telescopes. This kind of study has 
previously been performed in the low energy MSSM [|I], m. It is obvious that such 
a low energy approach offers more free parameters (in particular /i) to play with, 
and it thus becomes much easier to obtain heavy neutralino with a strong higgsino 
component. The main difference between the grand unification and the low energy 
approaches comes from the absence of neutralino annihilation into Higgs bosons for 
the GUT models, due to the RGE evolution of the SUSY parameters and to the mass 
spectrum hierarchy obtained. This possibility of lighter scalars (Higgs and squarks of 
first generation) in the low energy approach also sometimes enables fluxes from the 
Earth to be significant. However these low energy effective models are less consistent 
and do not take into account some nice features of SUSY GUT models such as the 
radiative electroweak symmetry breaking or the avoidance of Landau poles and CCB 
minima. 

We propose here to derive a low energy effective parameterisation of SUSY GUT 
models which appear to be favourable for the neutrino indirect detection of neutrali- 
nos. All these models present indeed the same characteristics: 

• a non- vanishing neutralino higgsino fraction J'h '■ 0.1 < ^ 0.4 (see figures ^ 
and I, I, 0), 

• XX W~^W~ , ZZ or XX ~^ tt dominant neutralino annihilation channels 
which are very efficient in order to get both a good value of the neutralino relic 
density and hard neutrino spectra. 

One notices that all the physics which governs these processes only involves neutrali- 
nos and/or charginos on top of standard particles. Therefore a minimal description 
assumes all scalars to be heavy and decoupled. The remaining necessary parame- 
ters are those defining the neutralino sector : Mi, M2, /i and tan/3. Moreover the 
unification relations and the RGE evolutions lead to M2 — 2Mi, while a neutralino 



higgsino fraction in the range 0.1 < /j;/ < 0.4 implies (see figure O) 



1.4 X lO-^M^ + 0.83Mi + 60 < < 1.5 X lO'^M^ + O.8M1 + 150 (5.1) 
with negligible dependence in tan/3 (see twin curves of figure |ll]). 



A SUSY GUT model satisfying at low energy: Ms^ai > Mi, M2 ~ 2Mi and n 
in the range indicated by equation ( |5ri| ) gets a cosmologically favoured value of the 
neutralino relic density and large muon fluxes coming from neutralino annihilation in 
the Sun. For a scan of low energy models with 50 < Mi < 1000 GeV, Mgcai ^ Mi, 
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Figure 11: Analytical iso higgsino fractions for tan/3 = 10 (lower twin curve) and 50 
(upper twin curve) shown as full lines, and the fitted functions described in equation ( IsTl] ) 
for fn = 0.4 (lower left) and //^ = 0.1 (upper right) shown as dashed lines. 



low energy parameterisation vs neutrino Indirect Detection 




Figure 12: Low energy parameterisation of RGE models versus experiment sensitivities. 
The thin full lines are the fitted functions of equation (5.2). 
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tan /3 = 5, 10, 20, 30, 40, 50, 60, and constrained by equation ( |5.1|) to get 0.1 < /// < 
0.4, we find the following range for the solar /i flux (see figure |1^) 

fH^0.1^5x 10^ — < log^oll^ Auxq] < 1.5 x 10^ — ^ ~ 0.4. (5.2) 

mi mi 

X. A. 

We thus see that independently of their high energy scale realisation, the most 
favourable models are within reach of the next generation (km^) of neutrino tele- 
scopes through the detection of neutrino fluxes from the Sun, up to masses of the 
order of 450 GeV for the most interesting relic densities. 



6. Conclusion 

In this paper, we have explored departures from CMSSM universality and singled out 
those most interesting for neutralino dark matter detection. The scalar sector can be 
of interest to adjust the neutralino relic density by opening sfermions coannihilation 
processes, but does not lead to increased detection rates because of the RGE evolution 
of first generation soft terms. The most determining parameters are the gaugino 
masses M2\gut and particularly M^Igut which respectively increase the wino and 
the higgsino content of the neutralino when lowered away from their universal values. 
The higgsino component is more efficient than the wino one to improve the detection 
rates, making of M^Igut the most relevant degree of freedom, as its value also affects 



the whole MSSM spectrum |^^. Such models with lower M^Igut values have a 
better relic abundance and are much more promising from a detection point of view, 
with rates increased by several orders of magnitude compared to the universal case. 
Some naturalness objections could be addressed on the high values of mo and mi/2 
commonly used in the literature. However it should be noticed that lower values of 



Ms tend to be more natural |H6[ . To finish, all the models with large neutrinos fiuxes 



studied here tend to have small SUSY contribution to the /x anomalous magnetic 
moment (a^^^"^ — 0), which might exclude them, but not before the value of the 
Standard Model itself is clearly ruled out. 
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